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Abstract. 

We study the influence of dipole-dipole interaction on the formation of vortices in 
a rotating dipolar Bose-Einstein condensate (BEC) of 52 Cr and 164 Dy atoms in quasi 
two-dimensional geometry. By numerically solving the corresponding time-dependent 
mean-field Gross-Pitaevskii equation, we show that the dipolar interaction enhances 
the number of vortices while a repulsive contact interaction increases the stability of 
the vortices. Further, an ordered vortex lattice of relatively large number of vortices 
is found in a strongly dipolar BEC. 
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1. Introduction 

The most fundamental properties of a superfluid is irrotational (V x v s = 0, v s the 
superfluid velocity) and resistance free flow. An important theoretical development 
came after the remarkable work by Landau that the excited state of a superfluid would 
represent the creation of vortices [lj. Quantized vortices have long been studied in 
superconductors, liquid Helium, etc., so called superfluids. The successful experimental 
realization of Bose-Einstein condensates (BECs) of alkali atoms in the past century has 
provided a new pathway to many developments in the understanding of the quantized 
vortices [2l[3llll[5l[6l[7l[8l[9j. In BEC experiments vortices have been nucleated with 
help of either by laser stirring or by rotating magnetic traps, and they are observed 
above a certain critical rotation frequency. The observation of quantized vortices could 
be considered as unquestionable evidence for the existence of superfluidity in BEC. 
Quantized vortices in BEC are having obvious analogy with that in liquid Helium and 
in type-II super conductors [TO] . 

Early experiments and theoretical studies on vortices in BEC have been mostly 
focused on the conventional Bose gas with local and isotropic interaction. However, 
many bosonic atoms and molecules have large dipole moments. For example, 
52 Cr, 164 Dy, and 168 Er BECs, which have a larger anisotropic long-range dipolar 
interaction superposed on the isotropic short-range atomic interaction, have recently 
been realized [HJ [121 EH EH US EU Ej . The study of dipolar BECs has revealed various 
interesting properties due to the peculiar competition between an isotropic, short-range 
contact interaction and an anisotropic, long-range dipolar interaction. These include, 
the dependence of stability on the trap geometry [HI CEj , new dispersion relations of 
elementary excitations [T8l [T9]. unusual equilibrium shapes, roton-maxon character of 
the excitation spectrum [2Ql EE EH [22] and novel quantum phases including supersolid 
and checkerboard phases [2H [25J [23] . 

Recently, there have been studies on vortices in dipolar BECs using mean-field 
models [23j. For instance, the structure and stability of vortices in dipolar BECs are 
found to be strongly affected by the anisotropic character of dipole-dipole interaction 
(DDI) and relative strengths of the dipolar and contact interactions [26j [271 [28] - 
Further, static hydrodynamic solutions, and dynamical stability and instability have 
been explored in a rotating dipolar BEC within Thomas- Fermi limit [29]. A second- 
order like phase transition of straight and helical vortex lines influenced by dipolar 
orientation has also been studied [30]. More recently the rotational properties of a 
dipolar BEC in a quasi two-dimensional (2D) geometry for an arbitrary orientation of 
the dipoles with respect to their plane of motion has been studied [3T] . 

However, most of the features of vortices in dipolar BECs have not been explored so 
far. Also there has been no studies on vortices in a strongly dipolar BEC. In the present 
paper, we analyze the formation of vortices in a dipolar BEC and the influence of strong 
dipolar interaction on the stability as well as on the number of vortices. By numerically 
solving the two-dimensional Gross-Pitaevskii (GP) equation, for a rotating trapped 
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dipolar BEC of 52 Cr and 164 Dy atoms in a quasi-2D geometry, we study formation and 
dynamics of vortices both in the presence and absence of s-wave (contact) interaction. 
The number of vortices in a pure dipolar BEC is found to increase when compared to 
that of a conventional BEC. We compare the vortex number estimated theoretically 
with that computed from numerical simulations. Further we notice that the inclusion 
of a repulsive contact interaction in a dipolar BEC enhances the formation of vortex 
lattices in the dipolar BEC. We also find that the critical rotation frequency decreases 
with the increase of the strength of dipolar interaction. 

The present paper is organized as follows. In Sec. [2] we provide a brief overview on 
the mean field Gross-Pitaevskii equation describing the properties of a rotating dipolar 
BEC confined in an axially symmetric harmonic trap potential. A two-dimensional 
reduction of the GP equation for a dipolar BEC under strong axial confinement is also 
discussed in Sec. [2j Then, in Sec. |3| we present the numerical studies on the formation 
and dynamics of vortices in dipolar BECs of 52 Cr and 164 Dy atoms. Here we analyze 
the formation of vortices in a pure dipolar BEC. We also explore the role of contact 
interaction on the stability, number of vortices and critical rotation frequency. Finally, 
in Sec. [4] we provide a summary and conclusion. 



2. Theoretical framework 



Many phenomenological properties of quantized vortices in BECs can be studied using 
mean field Gross-Pitaevskii (GP) equation [HI E] - A dipolar BEC with TV atoms, each 
of mass m at absolute zero temperature in a rotating frame can be described by the 
Gross-Pitaevskii equation as [26j 



. d<j>(r,t) 



iv 2 + V(r) + 4TvaN\(f)(r, t)\ 2 - QL Z 



+ N j Mr-Ol0(^)| 2 rfV]<KM), (1) 



where V(r) is the confining axially symmetric harmonic potential, 0(r,i) the wave 
function at time t with normalization J |0(r, t)\ 2 dr = 1, a the atomic scattering length, 
which can be tuned to a large extent via Feshbach resonance. The experimental value of 
the s-wave scattering length of 52 Cr atom is a = (103±4)a [33]. For the 164 Dy atoms the 
s-wave scattering length, yet to be measured, is assumed to be equal to that of 52 Cr [16J. 
The axial and radial trap frequencies are Q z cj and Vt p uj, respectively. In equation ([!]) 
length is measured in units of characteristic harmonic oscillator length / = ^h/muj, 
frequency in units of uo ) time t in units of uo~ l . L z = —i(xd y — yd x ) corresponds to the 
z-component of the angular momentum due to the rotation of the dipolar BEC about 
z axis with frequency Q. Here the Q is expressed in units of the radial trap frequency 
Q p cj. The integral term in equation ([!]) accounts for the dipole-dipole interaction with 

1 — 3 cos 2 9 

t/dd(x) = a dd rr^ , (2) 
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where x = r — r' determines the relative position of dipoles and 6 is the angle between 
x and the direction of polarization, z. The constant add = p^p 2 m / {12Tih 2 ) is a length 
characterizing the strength of dipolar interaction and, its experimental value for 52 Cr 
and 164 Dy are 15a and 130a , respectively [32], with a being the Bohr radius, ft 
corresponds to the magnetic dipole moment of a single atom and p the permeability of 
free space. 

The dimensionless three-dimensional (3D) harmonic trap is given by 

V(r) = \ny + \tf z z\ (3) 

where r = (p, z), with p the radial coordinate and z the axial coordinate. In a disk- 
shape, with strong axial trap frequency (Q z > Q p ), the dipolar BEC is assumed to be 
in the ground state, 

^ )= Hf^ eXP ("2^)' (4) 
of the axial trap so that the wave function 0(r) can be written as, 

0(r) = cj) 1D (z)(j)2D{x,y), (5) 

where 02d(^,?/) is the 2D wave function, d p = I and Xd 2 z = 1 with A = Q, z /Q, p the trap 
aspect ratio. Using ansatz ^ in equation ([T]), the z dependence can be integrated out 
to obtain the following effective 2D equation [341 1351 136] . 



dt 



+ V 2D - QL Z + -=-|0 2D (p,t)| 2 

72 



where V2D = [(1 + e)x 2 + (1 — e)y 2 ]/2 is the trap potential, e is a parameter used to 
introduce a small anisotropy for the nucleation of the vortices [33, [38] . In equation ([6]), 

n(K) = J exp(ik p .p)\(f)2D(p)\ 2 dp, (7) 

k p = (k x ,k y ), /&2d(£) — 2 — 3y / 7r^exp(^ 2 ) erfc(<^), and the dipolar term is written in 
Fourier space. 

In the following, we study the vortices of a disk-shaped dipolar BEC by solving 
the above equation ^ with the aid of a combined split step Crank-Nicolson and Fast 
Fourier Transform (FFT) based numerical scheme [35j[39]. The numerical simulations 
are carried out with dx = dy = 0.2 (space step) and dt = 0.005 (time step). 



3. Formation and dynamics of vortices in dipolar BEC 

Before studying the reduced GP equation (J6J), it is worth to comment on the validity 
of 2D approximation, at least, for the trap aspect ratios A = 30 and 100, which we are 
using throughout this paper. For this purpose, we compare the one-dimensional density 
profile of ([i]) and that from the numerical solution of full 3D equation ([T]) when Q = 0. 
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Figure 1. (Color online) One-dimensional density profile of a disk-shaped 52 Cr dipolar 
BEC, of 10 000 atoms from ID wave function Q, I^id^)! 2 and numerical solution of 
the 3D GP equation {l} with Vt = 0, |(/>^ m (>)| 2 , with trap aspect (a) A = 30 and (b) 
A = 100. The numerical solution of is obtained using imaginary time propagation 
with (a) dx = dy = dz = 0.2, dt 0.005 and (b) dx = dy = dz = 0.1, dt 0.0025. 



The numerical density profile is obtained by integrating out the x and y dependence of 
the full 3D density, 

/OO POO 
/ \<j>(x,y,z)\ 2 dxdy. (8) 
-oo J —oo 

In Figure [I] we plot the analytical and numerical density profiles of a disk-shaped 52 Cr 
dipolar BEC, of 10 000 atoms for A = 30 and 100. It is easy to see that the ID density 
profiles match very well. 

We prepare the ground state wave function by solving equation ^ numerically 
using imaginary time propagation in the absence of angular momentum = 0) and 
trap anisotropy (e = 0). The vortices are then created by evolving the ground state 
with the inclusion of angular momentum (Q ^ 0) in real time propagation up to 
t ~ 2 x 10 4 cc; _1 . A small anisotropy, e = 0.06 and a phenomenological dissipation 
(~ 10 -5 ) are introduced to facilitate the smooth vortex formation [371 HQ]. The 
presence of dissipative mechanism is evident in BEC experiments exhibiting collective 
damped oscillations [37] . Theoretical predictions based on mean-field GP equation with 
dissipative term found to give correct description of the damping. In a rotating BEC, 
vortex lattices will never be formed without dissipation even if the trapping potential is 
rotated fast enough. This phenomenological dissipation is effected by replacing the term 
Hd/dV with \i — rfd/dV in the time dependent equation ^ and 7 accounts for the 
dissipation. Vortices are formed for rotation frequencies above a critical value, ft > Q c . 
For smaller trap aspect ratios, a BEC with dominant dipolar interaction with the dipoles 
aligned along the trap symmetry axis becomes unstable [HI HT] . Here, we consider a 
dipolar BEC in disk-shape geometry, A 3> 1, as it is more stable. To explore the effects 
of dipolar interaction on the vortex formation we consider a disk-shape geometry with 
the number of atoms fixed to N = 10 000. 
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3.1. Vortices in a pure dipolar BEC 

First, we consider the case of a pure dipolar BEC with add = 15 a and a = with 
the trap aspect ratio, A = 30. The formation of vortex is analyzed by varying the 
rotation frequency Q in the range 0.1 to 0.9. When the condensate begins to rotate 
it elongates and the boundary surface becomes unstable. Then ripples are developed 
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Figure 2. (Color online) Contour plots of the density |</>2£>(/M)| 2 showing the 
development of vortices in a dipolar BEC with a = 0, a^d = 15 ao, N = 10 000, A = 30 
and ft = 0.6 after the trapping potential begins to rotate: (a) t = (b) t = 1000 cj _1 , 
(c) t = 2000 cj _1 , (d) t = 20 OOOcj -1 , (e) Phase distribution of condensate wave function 
</>2£>(/M) of fully developed vortices (steady vortex state) shown in (d), and (f) Three 
dimensional contour plot showing the vortices obtained from the numerical simulation 
of the full 3D equation ([I]) with the above parameters. 

on the surface and as a consequence vortices enter into the condensate above a critical 
rotation frequency, Q c ~ 0.38. As time progress these vortices approach to a stable 
configuration. The density plots of |02d(p,^)| 2 of the time development of vortices for 
Q = 0.6 are shown in Figures [2](a)-(d). About 7 vortices form a stable pattern as shown 
in Figure [2](d). We have calculated the phase from the final wave function £), which 
varies continuously from (dark) to 2tt (bright) and shown in Figure [2](e) . At the place 
where a vortex is located, a bifurcation (branching) in the phase is visible. 

We have also performed the numerical simulation using the full 3D equation ([!]) 
for the above set of parameters with the same rotation frequency, which again confirms 
the number of vortices and is shown in Figure [2](f). Further, the expectation value of 
angular momentum, defined as 

(L z ) =i / 4>*(p,t)(yd x - xd y )<j)(p,t) dp, (9) 
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Figure 3. (Color online) Time evolution of the expectation value of angular 
momentum (L z ) during the development of vortices shown in Figure [2] 



is also calculated as a function of time and is shown in Figure [3| The expectation value 
of angular momentum (L z ) gradually increases with periodic oscillation (see inset in 
Figure [3]) and then settles to almost steady value confirming the stable vortex pattern. 
The oscillation is due to the shape (quadrupole) deformation of the condensate at 
vortex free state and it is associated with creation of vortices. In Figures |4](a)-(c) 
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Figure 4. (Color online) Contour plots of the density distributions |</>2£>| 2 showing 
steady vortex state in a rotating dipolar BEC with a^d = 15 ao, a = 0, N = 10 000, 
A = 30: (a) Q = tt c = 0.38, (b) ft = 0.44, and (c) ft = 0.47. (d) - (f) The corresponding 
phase distributions of the condensate wave function <\>2D of (a) - (c). 
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and |ij[d)-(f) we depict the contour plots of the density |02d| 2 and the corresponding 
phase distribution of condensate wave function (j) 2 D for different rotation frequencies, 
Q = 0.38, 0.44, and 0.47, respectively. 

Next we consider the development of vortices in a dipolar BEC with larger dipolar 
strength, for example, add — 130 a . In this case, the critical rotation frequency is found 
to be Q c = 0.271, which is small when compared to that observed with add = 15 a 
above. Thus, a BEC with larger dipolar interaction nucleates vortices with smaller 
rotation frequency [28J. At critical rotation frequency, fi c a single vortex is formed as 
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Figure 5. (Color online) Contour plots of the density distributions |</>2£>| 2 showing 
steady vortices state in a rotating dipolar BEC with a^d = 130 ao, a = 0, N = 10 000, 
A = 30: (a) Q = tt c = 0.271, (b) ft = 0.29, and (c) tt = 0.33. (d) - (f) The 
corresponding phase distributions of the condensate wave function <p2D of (a) - (c). 



shown Figure [5](a), while the corresponding phase is depicted Figure [5](d). On increasing 
above Q c more number of vortices are created as illustrated in Figures [5](b)-(c) for 
Q = 0.29 and Q = 0.33, respectively, and the corresponding phases are shown in 
Figures [5](d)-(f). 

Further, it is interesting to note that in Figure [5| vortex arrays as well as a single 
vortex, are not centered with respect to the trap symmetry axis, in contrast to the 
centered vortex arrays that are shown in Figure [4] for a smaller dipolar interaction 
strength as it happens also in conventional BECs. In this case, the quadruple oscillation 
persists for longer times. This is essentially due to the strong dipolar strength, 
add = 130 ao, of 164 Dy atoms. 

On increasing Q further, the system with add = 15 a destabilizes as Q approaches 
unity. Moreover, a pure dipolar BEC with add — 130 a becomes unstable for Q > 0.62. 
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Rotation frequencies larger than 0.62 creates a large number of vortices and makes the 
boundary surfaces more unstable. Certain peculiar dynamics is observed during the 
formation of vortices when the rotation frequency Vt < 0.5 for larger dipolar strength 
(add — 130 a ). The dipolar BEC deforms as it rotates and the boundary surfaces 
become unstable, as a result more vortices enter and remains stable for a while. Then, 
as time progress, again the surfaces become unstable and few more vortices enter into 
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Figure 6. (Color online) Contour plots of the density |</>2£>(p, £)| 2 showing vortices in 
a rotating dipolar BEC with a = 0, a dd = 130 a , N = 10 000, A = 30 and tt = 0.35: 
(a) t = (b) t = 30 000 a;- 1 , (c) t = 40 000a;- 1 , (d) t = 50 000cj-\ (e) t = 80 000a;- 1 , 
and (f) t = 100 000 cj- 1 . 

the condensate and this scenario repeats for reasonably longer times. In this case (L z ) 
shows periodic oscillations for a very long time. In Figure [6] we plot the snapshots of 
density |02d| 2 showing the dynamics. Further, it is evident that the number of vortices 
increasing rather slowly as time progress. However, the above dynamics is not seen 
when the dipolar strength is relatively low or zero. 

We also examine the speed at which the vortices are formed in a dipolar BEC. 
In conventional BECs, with repulsive contact interaction, it takes about few hundred 
milliseconds for the propagation of surface ripples that are created due to trap 
deformation to the vortex formation [37l 138]. While in dipolar BECs, apart from trap 
deformation, the anisotropy nature of DDI breaks the axial symmetry more easily [3T] 
and stimulates rapid vortex formation. Our numerical simulations also confirm this 
rapid vortex formation as it take about few ten milliseconds for the creation of vortices, 
which is one order less compared to non-dipolar BECs. We considered a similar set of 
trap parameters, contact interaction strength and rotation frequency as in Ref. [37] . 
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Next, it is of interest to study the dependence of the vortex number (N v ) on the 
rotation frequency Q in dipolar BEC. In a rotating superfluid, the number of vortices 
formed for a given rotation frequency can be predicted using the Feynman's rule [42j, 
that is, N v = (mQ / ti) R 2 p (Q) , where R P (Q) is the radius of the rotating fluid. In the 
Thomas- Fermi (TF) limit, the radius of the trapped BEC under rotation is given by [7j, 

R P (Q) = R p (0) (l-- 2 ) . (10) 

Then the relation connecting vortex number and the rotation frequency can be written 

as 

N v = nR 2 p (o)(i-^J . (11) 

The Thomas-Fermi (3D disk/pancake) radius itW(O) is given by [22] 

R TF (0) = [15NX(a + 2a dd )} 1/5 . (12) 

Here we use R p (0) = R TF (0)/V3 @3]. 

The above relation (flO|) is used to study the vortex number in conventional BECs 



in a rotating deep optical lattice [44J. Here we use the above equation (11) to estimate 
the vortex number in a rotating dipolar BEC. This will be more useful to confirm the 
dominance of dipolar strength over Feynman's rule. We also calculate the vortex number 
as a function of the rotation frequency Q by numerically solving the the GP equation 
([6]). The numerical vortex numbers are calculated in the equilibrium state by evolving 
equation ^ in real time propagation for sufficiently longer times, typically of the order 
of 150 000CJ- 1 to 200 000 cj- 1 . 

In Figure [7|(a) we plot the vortex number N v obtained by using equation (11) 
(theory) and that calculated numerically (num) for a = 0, add — 15a and 130a for 
A = 30. We also show the vortex number for the case of conventional BEC with a^0 
and add — and dipolar BEC with a = 100a . For add — 15a the numerical and 
theoretical vortex numbers agree reasonably well. However, for larger dipolar strength, 
add — 130a , with higher rotation frequencies, Q > 0.6 the numerical results deviates 
from results of Thomas-Fermi approximation. The reason is due to the fact that for 
stronger dipolar interaction with larger rotation frequencies, the shape of the density 
profile may not be correctly described by the parabolic Thomas-Fermi profile. 

We also study the vortex number of a dipolar BEC confined in a harmonic trap 
with large trap aspect ratio, for example, A = 100. Figure [7|(b) shows the plots of 
numerical and theoretical vortex number N v with Q for a dipolar BECs with a = 0, 
add — 15a and 130a for A = 100. We find a similar dependence of the vortex number 
on the rotation frequency as above. However, the rotating dipolar BEC (add = 130 a ) 
is stable for larger rotation frequencies as against Q < 0.62 seen earlier with A = 30. 
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Figure 7. (Color online) Plot of the equilibrium vortex number (N v ) as a function 
of the rotation frequency Q for conventional BEC with a = lOOao, dipolar BECs with 
add — 15 ao ( 52 Cr) and a&d — 130 ao ( 164 Dy), and a = (pure dipolar) as well as 
a = lOOao for different trap aspect ratios: (a) A = 30 and (b) A = 100. Here 'theory' 
represents the vortex number estimated using equation ( pTj ) and 'num' corresponds to 
numerically computed values. 



3.2. Effect of contact interaction on vortices 

In the above, we have mostly analyzed the formation of vortices in a dipolar BEC 
with contact interaction set to zero. In general, the presence of a repulsive contact 
interaction enhances the stability of the dipolar BEC [14J. It is of interest to study the 
role of inclusion of contact interaction a on the formation of vortices. In order to see the 
effect of inclusion of contact interaction on the formation and stability of vortices we 
consider the case of rotating dipolar BEC with suitable repulsive contact interaction. 

First, we fix the contact interaction strength, say for example a = 10 a , and vary 
the Q to analyze the vortex formation in 52 Cr and 164 Dy condensates. In Figures [8](a)- 
(c) we plot the density profiles of the vortices with dipolar strengths add = 15 a ( 52 Cr), 
for different rotation frequencies Q = 0.5, 0.6 and 0.65, respectively. The corresponding 
phase distribution of the condensate wave function is plotted in Figures |8](d)-(f). For 
instance, a stable vortex lattice with eight vortices as shown in Figure |8](b) is found for 
add = 15 ao and Q = 0.6 with a = 10 ao. While in the absence of contact interaction 
it supports only 7 vortices as seen in Figure [2](e) . It is evident that the inclusion of 
repulsive contact interaction in a dipolar BEC increases the number of vortices. 

In Figures |9]( a) - (c) we show the contour plots of |02d| 2 of a stronger dipolar BEC 
with and a = 10a and add = 130a ( 164 Dy) for rotation frequencies Q = 0.5, 0.6 and 
0.65. Figures |9](d) - (f) depict the corresponding phase distributions of 2 d- As many 
as 18, 30 and 35, vortices for Q = 0.5, 0.6 and 0.65, respectively [Figures [9](a) - (c)] has 
been observed. We saw earlier that a dipolar BEC with a = and add — 130a becomes 
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Figure 8. (Color online) Contour plots of the density distributions \(/>2d\ 2 of a rotating 
dipolar BEC showing steady vortex state with a = 10ao, a^d = 15 ao, N = 10 000, 
A = 30: (a) Q = 0.5 (b) Q = 0.6 and (c) Q = 0.65. (d) - (f) The corresponding phase 
distributions of the condensate wave function <p 2 D of (a) - (c). 



unstable for rotation frequencies, Q > 0.62. Thus, as expected [H], the inclusion of 
a small repulsive contact interaction (a = 10a ) enhances the stability of the rotating 
dipolar BEC. 

Next we study the effect of dipolar strength in the formation of vortices for a 
fixed contact interaction strength and rotation frequency. By choosing the condensate 
parameters as a = 20 a , Q = 0.6, and A = 30 we analyze the pattern of vortices for 
different dipolar strengths. About 5 vortices were observed in the absence of dipolar 
interaction, that is, add — 0, as shown in Figure [l0]( a). This number increases to 10 



when add — 15a [Figure 10(b)] and further to 31 for add — 130a [Figure [10](c)] . This 
indicates that the dipolar interaction favors more number of vortices in a rotating BEC. 



We also calculate the critical rotation frequencies ^ c of a dipolar BEC with different 
contact and dipolar interaction strengths. The results are given in Table [T] The critical 
rotation frequency, for a rotating dipolar BEC decreases with the increase (or 
addition) of dipolar interaction or contact interaction. Similar conclusions have been 
drawn in earlier studies, which suggest that the critical rotation frequency depends on 
both the strengths of contact as well as dipole-dipole interactions [28]. However, the 
dipolar interaction play more dominant role in decreasing Q c than contact interaction. 
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Figure 9. (Color online) Contour plots of the density distributions |</>2£>| 2 of a rotating 
dipolar BEC showing steady vortex state with a = 10 ao, add = 130 ao, N = 10 000, 
A = 30: (a) Q = 0.5 (b) Q = 0.6 and (c) Q = 0.65. (d) - (f) The corresponding phase 
distributions of the condensate wave function <p 2 D of (a) - (c). 




Figure 10. (Color online) Contour plots of the density distributions \4>2d\ 2 of a 
rotating dipolar BEC showing steady vortex state with a = 20ao, ^ = 0.6, N = 10000, 
A = 30: (a) a d d = (b) a d d = 15 a and (c) a d d = 130 a . 



4. Summary and Conclusion 

We have studied the influence of dipole-dipole interaction on the formation of vortices 
in a rotating dipolar Bose-Einstein condensate of 52 Cr and 164 Dy atoms. In particular, 
we have performed numerical simulation of the quasi-2D Gross-Pitaevskii equation and 
showed that vortices appear more rapidly in a dipolar BEC when compared to non- 
dipolar BECs. Further, the number of vortices (N v ) is found to be significant in a 
dipolar BEC and increases monotonically with the increase of rotation frequency. We 
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Table 1. Critical rotation frequency Q c for vortices in a dipolar BEC 



a/ao 


add/ Q-o 







15 


0.380 





130 


0.271 


20 





0.500 


20 


15 


0.390 


20 


130 


0.275 


100 





0.360 


100 


15 


0.320 



also compared the number of vortices with Feynman rule. Our analysis suggests that for 
stronger dipole-dipole interaction strengths and higher rotation frequencies N v falls well 
outside the curve predicted by Feynman rule. We have calculated the critical rotation 
frequency (Q c ) for vortices, which decreases considerably with the increase of dipolar 
strength. 

The competition between the contact and dipolar interaction on the formation of 
vortices in a dipolar BEC is found to be complementary as the contact interaction 
enhances the stability of the vortices whereas the dipolar interaction increases the 
number of vortices. 
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